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Abstract — We consider a three-node half-duplex Gaussian 
relay network where two nodes (say a, b) want to communicate 
with each other and the third node acts as a relay for this two- 
way communication. Outer bounds and achievable rate regions 
for the possible rate pairs (R a , Rb) for two-way communication 
are investigated. The modes (transmit or receive) of the half- 
duplex nodes together specify the state of the network. A 
relaying protocol uses a specific sequence of states and a coding 
scheme for each state. In this paper, we first obtain an outer 
bound for the rate region of all achievable (R a , Rb) based on the 
half-duplex cut-set bound. This outer bound can be numerically 
computed by solving a linear program. It is proved that at any 
point on the boundary of the outer bound only four of the 
six states of the network are used. We then compare it with 
achievable rate regions of various known protocols. We consider 
two kinds of protocols: (1) protocols in which all messages 
transmitted in a state are decoded with the received signal in the 
same state, and (2) protocols where information received in one 
state can also be stored and used as side information to decode 
messages in future states. Various conclusions are drawn on 
the importance of using all states, use of side information, and 
the choice of processing at the relay. Then, two analytical outer 
bounds (as opposed to an optimization problem formulation) 
are derived. Using an analytical outer bound, we obtain the 
symmetric capacity within 0.5 bits for some channel conditions 
where the direct link between nodes a and b is weak. 

I. Introduction 

Two-way or bidirectional relaying has attracted significant 
interest recently Q] - © - In two-way relaying, a relay node 
assists in bidirectional communication between two nodes. 
Fig. [T] shows a three-node two-way relay network in which 




Fig. 1. Half-duplex Two-way Gaussian Relay Channel with Direct Link. 

nodes a and b want to communicate with each other at rates 
R a and Rb, respectively, while node r acts as the relay. While 
two-way relaying without the direct link between nodes a and 
b is studied in (lJ-Q, the more general two-way relaying 
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with the direct link is studied in [5|-|9|. In our work, we 
consider two-way relaying with the direct link. 

In (5}, achievable rate regions are derived for three proto- 
cols, namely the multiple access broadcast (MABC) protocol, 
time division broadcast (TDBC) protocol, and the hybrid 
broadcast (HBC) protocol under the decode-and-forward 
(DF) relaying scheme. In [6], the TDBC and MABC pro- 
tocols have been studied under other relaying schemes such 
as amplify-and-forward (AF), compress-and-forward (CF), 
mixed forward (MF) and Lattice forward (LF). A partial 
DF protocol, which is a superposition of DF and CF was 
studied in [10]. In |7), a three-phase cooperative MABC 
protocol (CoMABC), which outperforms the MABC and 
TDBC protocols in terms of sum rate in asymmetric channel 
conditions, was proposed. A transmission scheme based on 
doubly nested lattice codes was also proposed in JT), i.e., the 
CoMABC protocol uses the Lattice forward relaying strategy 
and is not a DF relaying protocol. 

The number of phases (or network states) in all the above 
protocols is between 2 and 4. Each phase (or network state) 
refers to a particular configuration of transmit and receive 
modes for the half-duplex nodes. In the absence of the 
direct link between a and b, the 2 states (multiple access 
and broadcast) in the MABC protocol are sufficient. In the 
presence of the direct link between a and 6, more states 
are required. The TDBC, HBC, and CoMABC protocols use 
more states to improve performance. However, all possible 
states have not been considered in these protocols. Since each 
node can be in transmit or receive mode, the 3-node half- 
duplex two-way relay channel has 2 3 (= 8) possible states. 
Of these 8, the two states in which all nodes are receiving 
or all nodes are transmitting are not useful in information 
transfer. Therefore, there are 6 useful states as shown in Fig. 
[2] A protocol that uses all these 6 states has been proposed 
in |8j, and the achievable rate region of this protocol has 
been derived. The improvement in achievable rate region with 
respect to the HBC protocol has been shown in |8). The 6- 
state protocol has also been presented independently in |9j 
in the context of a degrees of freedom analysis of two-way 
relaying with multiple-antenna nodes. 

In this paper, we first obtain an outer bound for the 
two-dimensional capacity region (of all possible rate pairs 
(R a , R b )) that is valid for all relaying protocols. This bound 
is derived based on the half-duplex cutset bound in JTT) and 
is more general than the bounds in |5|-[7|,which are for 
specific protocols - TDBC, MABC, HBC, and CoMABC. 



r 




Fig. 2. Possible states in a Half-duplex Two-way Relay Channel 

The outer bounds for TDBC, MABC, and CoMABC proto- 
cols consider at most 3 states and the bound for the HBC 
protocol considers 4 states. Further, the outer bound for the 
HBC protocol in J3J is difficult to compute and has not 
been computed. In contrast, the new outer bound can be 
numerically computed by solving a linear program. It is also 
proved that at any point on the boundary of the outer bound 
only four of the six states of the network are used. No outer 
bound is derived or compared with for the 6-state protocol 
in (8j. 

We then compare the new outer bound with achievable rate 
regions of the best known protocols so far - HBC, CoMABC, 
and 6-state protocols. We consider two kinds of protocols: 

(1) protocols in which all messages transmitted in a state 
are decoded with the received signal in the same state, and 

(2) protocols where information received in one state can 
also be stored and used as side information to decode mes- 
sages in future states. Through these comparisons, various 
conclusions are drawn on the importance of using all states, 
use of side information, and the choice of processing at the 
relay. In this context, a simple 6-state DF protocol without 
side information using all six states is also presented and 
compared with. 

Finally, two analytical outer bounds (as opposed to an 
optimization problem formulation) are derived. Using an 
analytical outer bound, we obtain the symmetric (R a = Rb) 
capacity within 0.5 bits for some channel conditions where 
the direct link between nodes a and b is weak. 

II. System Model 

Consider two nodes a and b communicating with each 
other in the network shown in Fig. [T] The relay node r 
assists this communication by receiving the information from 
these nodes and forwarding to the desired destination. All 
nodes are half-duplex nodes with a receiver noise variance 



of N. For simplicity, each node is assumed to have the same 
transmit power P for each state. Let h%, hi and /13 be the 
gains of channels a-r, b-r and a-b, respectively. The SNRs 

h 2 P h 2 P 

of these channels are denoted 71 = -j^-, 72 = -j^- and 

h 2 P 

73 = -h-. We use R a and Rb to denote the rate of data 
transmission (bits per channel use) from node a to node b and 
from node b to node a, respectively. We consider 73 < 71, 
73 < 72> i- e -' the direct link between a and b is weaker than 
the links a-r and b-r. Further, let 71 < 72, without loss of 
generality. Let C(j) = log 2 (l +7) represent the capacity of 
a complex gaussian channel with SNR of 7. 

III. Outer bound for any protocol 

In this section, we derive an outer bound for the capacity 
region of the half-duplex two-way Gaussian relay channel. 
This bound is an outer bound for any relaying protocol for the 
two-way relay network irrespective of the number of states 
and the relaying scheme (eg. AF, DF, CF) used at the relay 
node. This outer bound is derived using the half-duplex cut- 
set bound for information flow in an arbitrary half-duplex 
relay network in fTT) . In this two-way relay network, we 
have two flows, one from a to b and another from b to a. 
We derive the outer bound for the two flows considering two 
cuts that separate nodes a and b and information flow in both 
directions across the cuts. Combining these four bounds (2 
cuts x 2 directions), we get the rate region outer bound. 

Let Aj, i = 1, 2, . . . , 6, denote the fraction of channel uses 
in network state i. We obtain the outer bound in the following 
manner. (1) For a given real number k, upper bound the 
maximum possible rate R a subject to R a — kRb, and (2) 
Vary k and determine the whole region. For example, we 
can set k = tan# and vary 8 from to 90°. 

Theorem 1: Given R a = kRt, for some k > 0, the 
maximum possible Rb is upperbounded by Cbk obtained by 
solving the following linear program: 

C bk = max R b 

R b ,{\i} 

subject to 

kR b < AiC( 7 i + 73) + A 3 C( 7 i) + A 5 C (73) , 
kR b < \iC(rfa) + A 4 C( 72 ) + X 5 C ((V^+V^f) > 
R b < A 2 C( 72 + 73) + A 3 C( 72 ) + A 6 C (73) , 
Rb < A 2 C( 73 ) + A 4 C( 7l ) + A 6 C ((VtI + y/W) , 

6 

X] A < < 1. A * > 0. R b > 0. 

i=l 

Proof: For a half-duplex relay network with m states, for 
which the sequence of states is fixed with asymptotic fraction 
of time A; in state i, any achievable rate R of information 
flow is upper bounded as follows fTT| : 

ra 

R<mmJ2W(Xs;Y S c\X S c,i), (2) 



where a cut partitions the set of nodes into sets S and S c 
such that the source nodes are in S, the destination nodes are 
in S c , and S c is the complement of S. We use this bound 
for R a and R b as follows. In order to bound R a (rate of 
information flow from a to b), we consider the cuts defined 
by Si = {a}, and S 2 = {a, r}. Similarly to bound Rb, we 
use 5*3 = {b} and S4 = {b, r}. 

Using Si and S 2 as defined above, we get the following 
bounds. 

R a < mm{R a i,R a2 }, (3) 

where 

R al = A X I (X a ; Y r , Y b \i = 1) + A 3 / (X a ; Y r \X b , 1 = 3) 
+ X 5 I{X a ;Y b \X r ,i = 5), 



R a2 = XJ (X a ; Y b \i = 1) + XJ (X r ; Y b \i = 4) 



(4) 



(5) 



+ X 5 I(X a ,X r ;Y b \i = 5). 
Similarly, using S3, and S4, we get 

R b < mm{R b3 ,R b4 }, (6) 

where 

R b3 = X 2 I {X b ; Y r , Y a \i = 2) + A 3 / (X b ; Y r \X a , 1 = 3) 

+ X 6 I(X b ;Y a \X r ,i = 6) 



R M = X 2 I (X b ; Y a \i = 2) + XJ (X r ; Y a \i = 4) 



(7) 



(8) 



(9) 



+ X 6 I(X b ,X r ;Y a \i = 6) 

Now, we can upper bound the various mutual information 
terms as in fTT) to get: 

I{X a ;Y r ,Y b \i = l) <C( 7l + 73 ), 
I(X b ;Y r ,Y a \i = 2)<C( l2 +j3) ) 
I(X a ;Y b \X r ,i = 5)<C((l-p 2 5 ) l3 ), 
I(X b ;Y a \X r ,i = 6)<C((l-p 2 6 ) 73 ), 

where is the correlation coefficient between X a and X r 
and p 6 is the correlation coefficient between X b and X r . 

I (X a ; Y r \X b , 1 = 3) < C (71) , / (X b ; Y r \X a , 1 = 3) < C (72) 

(10) 

Similarily, we get 

7(X a ;y 6 |i=l)<C(73),/(X r ;Y 6 |i = 4)<C( 72 ) ! 

I{X b ;Y a \i = 2) <C( l3 ) 1 I(X r ;Y a \t = A) < C (71) , 
/ (X a , X r ; Y b \i = 5) < C (72 + 73 + 2p 5y /^) , 
I (X b , X r ;Y a \i = 6) < C ( 7l + 73 + 2p 6v /^) . 

(11) 

Using these mutual information upper bounds in the bounds 
^ and we get: 

R a < AiC( 7l + 73) + A 3 C( 7l ) + A 5 C ((1 - pi) 73) , 

R a < AiC(7 3 ) + A 4 C(7 2 ) + A 5 C (72 + 73 + 2p 5 %/7273~) , 

Rb < A 2 C( 72 + 73) + A 3 C( 72 ) + A 6 C ((1 - pi) 73) , 

Rb < A 2 C(7 3 ) + A 4 C(7i) + X 6 C (71 + 73 + 2 / o 6 \/7i73) , 

(12) 



Using C ((1 - p|) 73 ) < C (73) and C ((1 - pg)^) < C (73) 
in d9j, pg,/96 < 1 in ( fTT) , and i? a = fci?;,, we get the 
constraints in Q. Since these constraints have to be satisfied 
by any achievable rate for a given set of Aj's, the upper 
bound C a k can be determined by maximizing R a over the 
Aj's subject to the constraints of the linear program in ([T]). 

■ 

The complete outer bound for the capacity region can be 
obtained by solving the above linear program for different 
values of the ratio k. It is worth noting here that this outer 
bound is valid for any relaying scheme (since it is based only 
on the cut-set bound and mutual information bounds) for the 
two-way relay channel and is not restricted to DF or other 
protocols presented here. 

Using results from linear programming, we can now argue 
that for any point on the boundary of the outer bound (i.e., 
for any k), at most four of the six Aj's are non-zero, i.e., 
at most 4 states are used. However, it should be noted that 
which states are used depends on k. 

Theorem 2: For any k, there are at most four positive Aj's 
in the solution to the linear program in ([T}. 

Proof: It is optimal to choose Yli=i ^« = ^' ^ e -' satisfy 
the last constraint with equality. There are 4 other inequality 
constraints. Adding 4 slack variables, we can convert these 
constraints to equality constraints. Now, there are 7 + 4 = 11 
variables in the linear program with 5 equality constraints. 
Therefore, from | fT2l , we know that there are at most 5 
positive variables in the solution. Since R b should be positive, 
we can have at most 4 of the Aj's to be positive. ■ 

A. Alternative Outer Bound 

An alternative outer bound can be obtained as follows: 
(1) For a given real number k, upper bound the maximum 
possible weighted sum rate R a + kR b , and (2) Vary k and 
determine the whole region. 

Theorem 3: For some k > 0, the maximum possible 
R a + kR b is upperbounded by Ck obtained by solving the 
following linear program: 

Ck = max R a + kR b , subject to 

R a ,R b ,{\i} 

R a < AiC( 7 i + 73) + A 3 C( 7l ) + X 5 C (73) , 

Ra < AiC( 73 ) + A 4 C( 72 ) + A 5 C ((^72 + Vlsf) . 

Rb < A 2 C(7 2 + 73) + A 3 C(7 2 ) + A 6 C (73) , 

(13) 



Rb < A 2 C( 73 ) + A 4 C( 7l ) + A 6 C ((V7T + Vri) 2 ) 

6 

^Ai < 1, Aj > 0, R b > 0,;R a > 0. 



i=l 



For each k, we get a line. As k is varied, we obtain a set 
of lines that define an outer bound to the rate region. This 
outer bound is equivalent to the one obtained in Theorem 1 . 

IV. Comparisons 

Now, we compare the outer bound with the achievable rate 
regions of various protocols. We consider the following two 



classes of protocols: (1) Simple DF protocols with no side 
information across states, (2) protocols with side information 
across states. In addition to existing protocols, we also intro- 
duce a 6-state DF protocol with no side information across 
states and compare its achievable rate region with other 
protocols. This provides further insight on the importance 
of using all states and on the use of side information. We 
briefly describe the achievable rate regions here and then 
provide numerical results. 

A. DF protocols with no side information 

1) MABC protocol: The 3 -node half-duplex relay network 
has six useful states as shown in Fig. [2] The MABC protocol 
uses states 3 and 4. The achievable rate region for MABC, 
as been obtained in |5j, is described by the following 
inequalities: 



R a < min{A 3 C(7i), A 4 C(7 2 )}, 
R b < min{A 3 C(7 2 ), A 4 C(7i)}, 
Ra + Rb < A 3 C( 7 i+7 2 ). 



(14) 



2) 6-state DF Protocol: The proposed 6-state DF protocol 
uses all the 6 useful states. To optimize the fractions Ai, 
i = 1, 2, . . . , 6, we proceed as follows. Let Z l kl denote the 
rate of information flow (bits per channel use) from node 
k to node I in network state i. The flow rates Z l kl are 
constrained by the topology of the network state and the 
coding employed. States 3, 5 and 6 form Multiple Access 
(MAC) channels, while states 1 and 2 are Broadcast Chan- 
nels (BCs). For these states, we employ standard capacity- 
achieving coding methods fPJ) . In state 4, the relay node 
r broadcasts the information it received from nodes a and 
b in previous states. Thus, in state 4, each receiver already 
knows the message that it transmitted to the relay, i.e., state 
4 is a BC with two receivers knowing the message intended 
for the other user a priori. The capacity region for state 4 
was determined in fT4) , and we employ the coding scheme 
given there. 

The achievable rate region of the half-duplex bidirectional 
relay channel with the 6-state DF protocol is the closure 
of the set of all points (R a ,R b ) satisfying the following 
constraints: 



Ra - Z a 



z. 



ab 



7 2 



7 2 
^ba 



' ^ab 
Zba 



' 3 

ar' 

3 

br • 



^Ai = 1,0 < A, < 1,0 < at,a 2 < 1 

State 1 rate constraints: 

Zl r <\rC{a^i),Z\ h <\ x C 
State 2 rate constraints: 

Zlr < A 2 C(0!272), Zl a < \ 2 C 

State 3 rate constraints: 



(1 - ai)7 3 
1 + ai7 3 

(1 - Q!2)73 
1 + «273 



Z'l < A 3 C( 7 i), Zt < A 3 C( 72 ), Z'l+Zl < A 3 C ( 7 i + 72) 



State 4 rate constraints: 

Z A ra < A 4 C( 7 i), Z^ b < A 4 C( 72 ). 

State 5 rate constraints: 
Zl h < A 5 C( 72 ), Z 5 ab < \ 5 C( l3 ),Z 5 rb +Z 5 ab < A 5 C (72 + 73) ■ 

State 6 rate constraints: 

Z e ra < A 6 C( 7l ), Zt < A 6 C( 73 ), Z e ra +Z 6 ba < X 6 C (71 + 73) ■ 

Flow Constraints: Information received at node r from 
node a should be equal to the information forwarded from 
node r to node b. Similarly in the other direction also. Thus, 
we get the following equality constraints: 

Zlr + %ar = ^rb + ^rbi %br + ^6r = ^ra + ^rer 

In these constraints, a± is the fraction of power used for 
the message from a to r in state 1, and a 2 is the fraction 
of power used for the message from a to r in state 2. 
Numerical evaluation of this achievable rate region of the 
6-state DF protocol is done for some illustrative examples in 
the numerical results section to show the importance of the 
various states in achieving different parts of the rate region. 

B. Protocols with side information across states 

1) TDBC and HBC protocols: TDBC [5] is a three phase 
protocol in which states 1, 2 and 4 are used. In first phase 
(state 1), a transmits at a rate equal to the capacity of 
link between a and r. At this time node b listens to this 
transmission and uses this information as side information 
for decoding after phase 3 (state 4). Phase 2 is similar to 
phase 1 in which b transmits at a rate equal to the capacity 
of link between b and r and a listens to this. At the end of 
each of these phases, relay node r decode the messages and 
does a binning operation on these messages. In phase 3, r 
transmits this binned information to a and b. Since a and 
b know the message meant for the other destination, relay 
node r can use XOR of these messages for broadcasting. 

HBC is a four phase protocol in which states 1, 2, 3 and 
4 are used. States 1, 2 and 4 are used in the same way 
as in TDBC. State 3, where terminals a and b transmit 
simultaneously (MAC) to relay r, is also added. These 
messages are decoded at the relay and forwarded as such 
in state 4. In HBC, state 4 is used for forwarding binned 
messages from state 1 and 2 as well as for forwarding 
messages received in state 3. The HBC protocol is always 
better than TDBC and MABC protocols since they are special 
cases of the HBC protocol. The use of side information 
from one phase in decoding during another phase provides 
improvement in TDBC and HBC over MABC for some 
channel conditions. 

The achievable rate region for HBC protocol has been 
obtained in |5|. It is the closure of the set of all points 



(R a ,Rb) satisfying following constraints. 

R a < min{(A 1 + A 3 )C( 7 i), AiC( 73 ) + A 4 C( 72 )}, 
R b < min{(A 2 + A 3 )C( 72 ), A 2 C( 73 ) + A 4 C( 7l )}, 
R a + Rb< AiC( 7 i) + A 2 C( 72 ) + A 3 C( 7 i + 72 ), (15) 



4 

£ 



A; 



1,0 < \ < 1. 



2) 6-state protocol: States 1, 2, 3 and 4 are used in the 
HBC protocol. A 6-state protocol that uses all 6 states has 
been proposed in (8). This protocol is similar to the HBC 
protocol except that MAC states 5 and 6 are used before the 
use of state 4. Here, we arrive at the same achievable region 
as that of the 6-state protocol in a slightly different way than 
in (SJ. States 1, 2, 3 and 4 are used as in the HBC protocol. 
In states 5 and 6, we have MAC transmissions from {a,r} 
and {b,r} to node b and a, respectively, i.e., in states 5 and 
6, there is a direct transmission between the terminal nodes 
a and b, and forwarding of messages received at relay in 
previous phases. In the HBC protocol, the direct link between 
a and b is only used to obtain side information to decode the 
transmission from the relay. In the 6-state protocol, the direct 
link is also used to send messages. 

The achievable rate region of the 6-state protocol is now 
the closure of the set of all points (R a ,R b ) satisfying the 
following constraints: 

i?a<^fc + min{(Ai+A 3 )C( 7l ), 

AiC( 73 ) + A 4 C( 72 ) + Z?, b }, 
min{(A 2 + A 3 )C( 72 ), 



Rb < Z\ 



ba 



A 2 C( 73 ) + A 4 C( 7l ) + Z 6 ra }, 



R a + Rb< AiC( 7l ) + Zl b +\ 2 C( l2 ) + Zt a + A 3 C( 7l + 72 ), 

A,: = 1,0 < K < 1. 
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sources. Thus, in state 3, it is bounded by the MAC capacity 
bound, and in other states it is bounded by individual links. 
In state 4, there is no information flow from source. 

Now, we can observe that the choice of Z^ b , Z^ b , Z ba , 
and Z® a that maximize the achievable rate region are: Z^ b + 
Zl h = A 2 C( 72 + 73 ), Zt + Z 6 ra = A 6 C( 7l + 73 ), Z\ h = 
A5C( 73 ), and Z ba = A 6 C( 73 ). Substituting this choice, we 
get the same region as in (8). 

3) CoMABC protocol: The achievable rate region for 
CoMABC protocol is taken from ITJ. We are assuming 
7 3 < 7i < 72- Thus, states 3, 4 and 6 are used in this 
protocol. The constraints for the CoMABC protocol are given 
by 

R a < min{A 3 i?* r , A 4 C( 72 )}, 

R b < mm{X 3 R* br + A 6 C( 73 ), A 4 C( 7i ) + A 6 C( 7 i + 73 )}, 
7i 



Z 5 ab < A 5 C( 73 ), Z 5 rb < A 5 C( 72 ), Z 5 ab + Z 5 rb < A 5 C( 72 + 73 ). 
Z 6 ba < A 6 C( 73 ), Z 6 ra < A 6 C( 7l ), Zt + Z e ra < A 6 C( 7l + 73 ). 

The above achievable rate region is obtained by modifying 
the achievable rate region for HBC protocol to include 
the effect of states 5 and 6 as well. The main points are 
summarized here. (1) Z\ b represents the direct information 
transmission from node a to node b in state 5. Similarly Z ba 
represents the direct information transmission from node b 
to node a in state 6. (2) The first term inside min{} in 
the constraints for R a and Rb correspond to the flow from 
source to relay and the second term in the constraints for R a 
and R b correspond to the flow from relay to destination. For 
source to relay flow, this protocol also uses the same states as 
HBC protocol. For forwarding from relay, this protocol uses 
states 5 and 6 also. Thus, while the first terms inside min{} 
remain the same as HBC, the second term has an additional 
term corresponding to flow from relay in states 5 and 6 for 
rates R a and R b , respectively. (3) The bound on R a + R b 
corresponds to the sum of all information flow from both 



7~>* 



R br — 



log 
log 



7i + 72 

72 

,7i + 72 
C. Numerical Results 
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Fig. 3. Achievable rate region and Outer bound for various protocols: 
71 = 10 dB, 72 = 15 dB, 73 = 3 dB 

1) Comparison with DF protocols: Figs. [3]and|4]show the 
comparison of different DF protocols with the outer bound. 
The 6-state DF and HBC protocols are always better than 
MABC. The 6-state DF protocol also achieves several rate 
pairs that the HBC protocol cannot achieve even though the 
HBC protocol uses side information across states. The HBC 
protocol also achieves some rate pairs that the 6-state DF 
protocol cannot achieve (see Fig. [4] for small R b ). The 6- 
state protocol achieves a larger rate region than all other 
protocols as expected since it uses all states as well as the 
side information used in HBC protocol. The achievable rate 
region of the 6-state protocol is closer to the outer bound in 
Fig. [3] where the SNRs are higher than in Fig. |4] 

2) Comparison with CoMABC protocol: Figs. [5] [6] and [7] 
show the comparison of the 6-state DF and 6-state protocols 



□ 


6-state DF 




MABC 


< 


HBC 





6-state 


A 


Outer Bound 





Fig. 4. Achievable rate region and Outer bound for various protocols: 
71 = dB, 72 = 5 dB, 73 = -7 dB 



Fig. 6. Comparison with CoMABC: 71 = 30 dB, 72 = 35 dB, 73 = 13 
dB 



with the CoMABC protocol. As mentioned earlier, CoMABC 
is not a DF protocol and the relay forwards a estimated 
function of the two messages. Thus, it has some advantage 
over DF protocols. However, the CoMABC protocol does 
not use all states. Therefore, it does not achieve some rate 
pairs that 6-state or 6-state DF can achieve. The CoMABC 
protocol is optimized to maximize sum rate. Therefore, it 
performs well near the maximum sum rate points. States 5 
and 6 used in the 6-state DF and 6-state protocol are very 
useful for achieving asymmetric rates. Fig. [6] considers a high 
SNR scenario where 71 = 30 dB, 72 = 35 dB, and 73 = 13 
dB. In this case, the CoMABC protocol almost achieves the 
outer bound for symmetric rates, while the 6-state protocol 
is close to the outer bound for asymmetric rates. Fig. [7] 
considers a low SNR scenario where 71 = dB, 72 = 5 
dB, and 73 = —7 dB. In this case, the CoMABC protocol is 
worse than the 6-state protocol except near R a = 0. 




Fig. 5. 
dB 



Comparison with CoMABC: 71 = 10 dB, 72 = 15 dB, 73 



V. Analytical Outer Bounds 



The outer bound in Theorem [T] can be evaluated by solving 
a linear program for each k. In this section, we derive an 
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Comparison with CoMABC: 71 = dB, 72 = 5 dB, 73 



analytical expression that is an outer bound by analyzing the 
dual of the linear programs in Theorems [T] and [3] These outer 
bounds are easier to evaluate and can also be used to gain a 
better understanding of the bound. 
A. Outer Bound 



The dual of the linear program in Theorem |Tj 

subject to y 5 > 2/iC(7i + 73) + ViC (73) , 
2/5 > 2/3^(72 + 73) + 2/4C (73) , 
2/5 > yiC(7i) +y 3 C (72) , 
2/5 > 2/2C(7 2 ) + 2/4C (71) , 
2/5 > 2/1^(73) + 2/2C ((V72 
2/5 > 2/3C(73) +2/4C ((VTi 
kyi + ky 2 + 2/3 + 2/4 > 1. 
Theorem 4: For any k > 1, Rb is upperbounded as 

R b < max{Ti,r 2 ,r3,T 4 }, 
3fc-l C( 7 i)C( 72 ) 



is: mini/5 

{Vi} 



(17) 



73 J 
73? 



where T\ = 



2fc 2 C( 7l )+C(7 2 ) ! 



To = 



2k - 1 C( 72 )C( 7 i + 73) + C( 7 i)C( 73 ) 



To 



T 4 



2fc 2 C( 7l )+C( 72 ) 

2fc - 1 C( 72 )C( 73 ) + C( 7 i)C((V7l + V73O 2 ) 
2fc 2 C( 7l )+C( 72 ) 

1 C( 7 i)C( 73 )+C( 72 )C((V 7 T+y 7 i) 2 ) 
2fc C ( 7l )+C( 72 ) 



Proof: The value of the dual program at any feasible 
point is an upper bound on the value of the primal problem. 
Choosing an appropriate feasible point in the dual program 
provides a good upper bound on Rb- The bound on Rb is 
obtained by choosing the following: 

2 k - 1 C( 72 ) = 2fc-l C( 7l ) 

C( 7l )+C( 72 )' 2/2 2fc2 C( 7l )+C( 72 )' 

_ 1 C( 7 i) . _ 1 C( 72 ) 



2/1 



2/3 



2fcC( 7l )+C( 72 )' 



.'74 



2fcC( 7l )+C( 72 ) 



Corollary 1: (1) For fc = 1, T 2 < T 4 . Therefore, R b < 
max{Ti, T 3 , T4}. (2) For 7 i = 72 = 7 and k = 1, we get 

< max I ^M, 1 [C( 73 ) + C((V7 + VtsT)] } ■ (18) 

Remark 1: A result similar to Theorem [4] can be obtained 
for k < 1 as well. In this case, we can set Rb = k'R a , where 
k' > 1 and use the same technique as in Theorem |4] The 
expressions obtained for T\ to T4 are similar to Theorem [4] 
except that 7 i and 72 are interchanged in each expression. 

Remark 2: Rb is also upper bounded by the upper bound 
for one-way relaying from & to a with R a = 0, i.e., we have 



chi + 73 )C((V72~ + VtsT) - c2 (73) 



C( 7 i + 73 ) + C((V7S + v/Ts) 2 ) - 2C (7s) 



R h < 



This bound is obtained by solving the dual program for one- 
way relaying. When there is no direct link, i.e., 7 3 = 0, this 
reduces to Rb < ^(j^+cjyl) ' wriere me bound is the capacity 
of the half-duplex two-hop linear network pT) . 

B. Capacity Results 

Theorem 5: (1) For 71 = 72 = 7 and k = 1, the upper 
bound on R a is C( 7 )/2 for 73 < 73 o, where 73 o satisfies 

/( 730 ) = 2C( 7 ), /( 73 ) = C( 73 ) + C((V7 + ^7i) 2 ). 
(2) For 7 i ^ 72 and k — 1, the upper bound on i? a is 
C( 7 i)C( 72 )/(C( 7 i) +C( 72 )) for 73 < min( 7 3i, 732 ), where 
731. 732 satisfy /i( 73 i) = / 2 ( 732 ) = 2C( 7l )C( 72 ) where 

/i( 73 ) = C( 72 )C( 73 )+C( 7l )C((V 72 -+ N /7i) 2 ) and / 2 ( 73 ) = 
C(7i)C( 73 ) +C( 72 )C((y 7 T+ y^) 2 ). 

Proof: Note that /(•), /i(-), and / 2 (-) are monotonically 
increasing functions of 73 for a given 7 i, 72 . Therefore: (a) 
unique solutions exist for 73 o, 73 i, and 732 , and (b) the 
resulting max term in each upper bound in Corollary [T] is 
the first term. ■ 
In the 7l = 72 = 7 , k — 1 case, the above result 
means that the upper bound is C( 7 )/2. The lattice coding 
scheme in [151 can achieve rates within 0.5 bits of this 



upper bound without any direct link. Therefore, for the weak 
direct link regime specified by 73 o, C( 7 )/2 is the capacity 
within 0.5 bits and can be achieved without using the direct 
link. In the 7l 7^ 72 , k = 1 case, the bound reduces to 
C( 7 i)C( 72 )/(C( 7 i) +C( 72 )) when the direct link is weak. 
This is equal to the one-way relaying bound when 7 j 7^ 72 
without the direct link. 

C. Analytical Upper Bound for Weighted Sum Rate 
Theorem 6: For any k , kR a + Rb is upperbounded as 

kR a + Rb < max{Ti, T 2 , T 3 , T 4 }, 

rp C( 7 i+73)C((V72 + V7^) 2 )-C 2 (73) 

W6ie 1 C( 7l+73 )+C(( % / 72 -+V7i) 2 )-2C( 73 )' 

C( 72 + 7 3 )C((VtT+ VTsT) - c2 (7s) 



T 2 = 



T, = k 



+ 



C( 7 2 + 73) + C((VtT + V73~) 2 ) - 2C( 73 ) ' 
C( 7 i)C(V72"+/7i) 2 )-C(7 3 ) 



C(ti + 73) + C(V72 + V73) 2 ) - 2C( 73 ) 
C( 72 )C(V 7 T+ J^) 2 ) -C( 73 ) 



T 4 = fc 



C( 72 + 73) + C(VtT + V73) 2 ) - 2C( 73 ) ' 
C(7 2 )C(7i +73) -C( 73 ) 



C(ti + 73) + C(V72" + VTs) 2 ) - 2C( 73 ) 

^(71)^(72 +73) ~ C( 73 ) 

C( 7 2 + 73) + C( v AyT + 7^)2) - 2C( 73 ) ' 



Proof: The bound is obtained by writing the dual of the 
linear program and choosing the dual variables as: 
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2/2 = 



2/3 
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C(7i + 


73) - 
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l-VTi) 2 )- 
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C(72 + 
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C( 72 - 


f73)+C(( A 


/tT- 


t-VTi) 2 )- 


2C( 73 ) 



D. Numerical Results 

Figure [8] compares the numerical (solving a linear pro- 
gram) and analytical outer bounds for three channel con- 
ditions A, B, and C. The two analytical bounds are close 
to the numerical bound. The first analytical bound is better 
for symmetric rates, while the second outer bound is better 
for asymmetric rates. In case B, where 7l = 72 , the first 
analytical bound matches with numerical bound for k = 1 
as expected from Theorem [5] Figure [9] shows the threshold 
min( 73 i, 732 ) defined in Theorem [5] as a function of 72 for 
a given c = 7 i/ 72 . Three values of c are considered. The 
results illustrate for these cases the range of 73 < min( 73 i, 
732 ) when the direct link can be ignored in the outer bound. 
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Fig. 8. Comparison of analytical and numerical outer bounds - Case A: 
7i = 10 dB, 72 = 15 dB, 73 = 3 dB, Case B: 71 = 20 dB, 72 = 20 dB, 
73 = 8 dB, Case C: 71 = 30 dB, 72 = 35 dB, 73 = 13 dB, 
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Fig. 9. Threshold for 73 = min(73i, 732) with 71 = C72 



VI. Summary 

The half-duplex two-way Gaussian relay channel with 
direct link has been studied. First, an outer bound for the 
rate region achieved by any protocol was derived. At any 
point on this outer bound, only four out of the six possible 
states are required. Then, this outerbound was compared with 
various existing protocols - HBC, 6-state, and CoMABC 
- and their achievable rate regions. While the CoMABC 
protocol is good near the maximum sum rate region, the 
6-state protocol is good for asymmetric rates, particularly 
at lower SNRs. The 6-state protocol is able to achieve this 
even while being restricted to only DF. The analytical outer 
bounds derived in this paper were shown to be close to the 
numerically computed outer bound. These analytical outer 
bounds can be used to gain a better understanding of the 
capacity region boundary. We also obtained the symmetric 
capacity (to within 0.5 bits) using the analytical bounds for 
some channel conditions where the direct link between nodes 
a and b is weak. The analytical bounds using the dual of the 
linear program could be potentially improved in the future 



to identify more channel conditions where capacity can be 
determined. 
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